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摘 要 : 为 精确 识别 雪崩 流动 特性 及 流 态 信息 ,全 


面 分 析 其 运动 过 程 。 本 研究 基于 无 人 机 倾斜 摄影 技术 获取 高 分 


状 率 航拍 数据 ,以 阿尔 先 沟 雪崩 易 发 区 为 例 ,通过 现场 调查 .无 人 机 遥感 解 译 精细 探测 雪崩 活 动 过 程 ,确定 RAMMS 


模型 输入 参数 ,在 此 基础 上 对 不 同类 型 雪崩 事件 进行 模拟 和 重建 ,对 比分 析 传 统 地 
模拟 结果 的 差异 ,探讨 不 同类 型 不同 雪 层 释放 条 件 下 雪崩 活动 过 程 。 人 研究 结果 表 
而 调查 方法 结合 无 人 机 遥感 和 数值 模拟 相互 验证 ,提高 了 灾害 发 育 状 况 评 佑 的 
准确 性 。(2) 2 月 中 名 阿尔 先 沟 坡 面积 雪 厚 度 趋 近 于 临界 厚度 值 ,持续 降雪 使 雪 层 失 稳 触 发 新 雪 雪 前 。 调 查 时 仍 处 
下 雪 榴 自重 逐步 增 大 ,有 超过 雪 的 抗 断 强度 的 趋势 ,整体 稳定 性 较 
,释放 量 可 达 8.2669x10' mi’, 运动 时 长 约 为 128 s, 并 在 120 
这 动 高 度 达到 最 大 , 约 为 3.55 m, 最 大 流动 速度 为 18.34 ms', 最 大 冲击 力 可 达到 32.67 kPa, 形 成 面积 
积 体 。 通 过 相互 验证 , 坡 


的 雪崩 调 查分 析 体 系 ,将 传统 的 地 


阶段 , 雪 层 裂缝 加 剧变 形 , 风 力作 / 
量 上 方 积 雪 平台 为 潜在 释放 区 的 坡 面 型 雪崩 


于 灾害 孕育 
差 ,(3) 以 坡 
s A HERA 
3369.7 m°, (RFR 1.8525x10" m° hY 


下 调查 IC) LR EARS 
明 :(1) 以 倾斜 摄影 技术 为 核心 


看 型 雪 衣 并 非 积 雪 平 台 的 释放 ,地 面 调查 结果 与 数值 模 


拟 
右 


解 译 结果 存在 差异 (4) 沟 槽 一 坡 二 
, 雪 骨 持续 时 间接 近 于 300 s, 


取 的 准确 性 ,可 为 今后 雪崩 潜在 危险 预测 风险 规避 及 灾害 
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雪崩 是 气象 . 雪 特 性 、 地 质地 貌 等 因子 协同 作 
用 下 ,山坡 积 雪 失 稳 引起 的 严重 自然 灾害 ,其 不 仅 
阻碍 着 区 域 交 通 .旅游 业 .农林 牧 生产 发 展 ,而 且 威 
胁 着 人 民生 命 财产 安全 ,给 社会 造成 了 巨大 压力 ? 。 
随 着 山区 经 济 的 快速 发 展 ,我 国 对 雪崩 防治 工作 提 
出 了 更 高 要 求 “”。 而 雪崩 调查 研究 工作 是 关系 防 
治 工 作 科学 性 和 有 效 性 的 基础 ,因此 蝇 待 开展 雪崩 
调查 人 研究 ,以 更 详细 和 准确 的 调查 资料 支撑 决策 者 
对 于 雪崩 风险 评判 。 

雪 骨 运动 轨迹 运动 特征 等 是 灾害 防治 的 关键 
因素 。 然 而 ,目前 针对 大 规模 雪崩 的 调查 研究 , 受 
限于 高 寒山 区 恶劣 自然 条 件 地形 条 件 和 交通 可 达 


收 稿 日 期 : 2023-05-19; ”修订 日 期 : 2023-07-11 
基金 项 目 : 交通 运输 行业 重点 科技 项 目 (2022-ZD6-090); 新 疆 交通 运输 科技 项 


i 复合 型 雪崩 为 沟 权 坡 
E 积 区 最 大 流动 速度 6.58 ms ,最 大 冲击 力 17.97 kPa, 平 均 堆积 深度 为 1.64 m, 影 响 
围 1178.5m2, 堆 积 量 3107.76 mi ,地面 调查 结果 与 数值 模 ] 


折 雪 层 断 裂 释 放 且 断裂 深度 仅 为 临界 厚度 值 的 60% 左 


以 结果 一 致 。 人 研究 结果 一 定 程度 上 提高 了 雪崩 事 件 信 息 获 
应 急 处 置 提供 强 有 力 的 数据 支撑 及 科学 依据 。 


性 的 影响 ,传统 的 人 工地 面 调查 方式 往往 不 能 
面 . 客 观 地 对 雪崩 运动 轨迹 做 出 准确 描述 且 会 耗费 
大 量 的 人 力 和 时 间 成 本 ”“ 。 因 此 ,迫切 需要 发 展 能 
够 快速 .全面 精准 描述 雪崩 运 动 过 程 的 新 技术 。 
近年 来 ,无 人 机 航 摄 系统 因 其 高 时 效 、 高 分 辨 
率 、 低 成 本 、 低 损耗 \ 低 风险 及 可 重复 等 诸多 优势 ， 
被 越 来 越 广 泛 地 应 用 于 灾害 调查 等 领域 ”。 和 鲁 恒 
等 "提出 无 人 机 获取 的 影像 在 三 维 可 视 化 中 能 提供 
更 详细 ,更 丰富 的 几何 和 语义 信息 。 赵 星 涛 等 "在 
无 人 机 低空 摄影 技术 获取 的 高 分 辨 率 正 射影 像 及 
大 场景 立体 正 射 影像 基础 上 对 矿山 地 质 灾 害 进 行 
建 模 分 析 ,精细 探测 灾害 活动 过 程 并 定量 与 定性 分 
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析 矿 山地 质 灾害 的 类 型 .规模 、 损 失 及 危害 等 。 芯 
大 雷 等 "在 利用 无 人 机 低空 摄影 技术 中 发 现 其 能 
很 好 地 认识 灾害 空间 分 布 规律 和 发 育 特征 , 且 能 清 
楚 认 识 灾害 变形 迹象 和 成 灾 过 程 。Eirik 等 "将 无 人 
机 航拍 获取 的 雪崩 运动 轨迹 正 射 影像 作为 验证 数 
据 , 在 一 定 程度 上 提高 了 评估 高 分 辩 紊 的 SAR 用 于 
探测 雪 朋 运动 的 准确 度 。Yang 等 ”将 无 人 机 作为 技 
术 手 段 发 现 航 拍 数据 能 够 更 好 的 识别 雪 骨 形态 特 
征 。Markus 等 “学 者 利用 无 人 机 进行 了 雪 朋 空 
调查 ,根据 获取 的 雪 毅 碎片 正 射 照 片 ,提取 到 了 有 
关 流动 特征 和 流 态 的 信息 ,并 生成 雪 骨 表面 DEM 数 
据 ; 该 数据 能 够 更 有 效 的 运用 于 雪崩 动力 学 建 模 和 
雪 深 空间 变异 性 研究 。 相 较 于 传统 人 工 实地 调查 
方式 ,上 述 非 接触 式 获取 灾害 信息 的 技术 在 精度 、 
效率 安全 性 等 方面 都 有 着 明显 的 优势 ,可 为 保障 
人 民生 命 财 产 安 全 和 灾害 应 急 响 应 等 工作 提供 更 
充分 的 科学 依据 。 

在 雪 朋 运动 过 程 模拟 领域 ,由 于 雪 朋 的 突 发 性 
及 运动 过 程 复 杂 , 现 有 的 模拟 模型 大 多 基于 一 些 经 
验 模 型 及 简化 的 物理 模型 “。 此 类 模型 主要 针对 
运动 过 程 中 的 特征 值 进行 计算 。 而 瑞士 联邦 冰雪 
与 雪崩 研究 所 在 Voellmy" fI Salm 等 "所 提出 的 
Voellmy-Salm( 简 称 VS) 模 型 的 理论 基础 上 研发 了 二 
维 雪 有 衣 动 态 数 值 模 型 
Simulation (简称 RAMMS)。 该 模型 考虑 雪 骨 运动 过 
程 中 流动 高 度 和 速度 变化 的 非 定 值 和 非 均匀 运动 
特征 ,使 用 积 雪 深 度 平 均 方程 同时 耦合 随机 动能 
(Random Kinetic Energy , RKE ) 模 型 来 考虑 质量 和 动 
能 守恒 ,解释 了 雪 粒 之 间 的 随机 运动 和 非 弹 性 互相 
作用 ,从 而 能 够 更 好 的 模拟 运动 过 程 中 各 特征 参数 
的 动态 变化 趋势 。 因 此 ,为 了 提高 模型 计算 结果 的 
准确 性 ,输入 参数 将 更 依赖 于 高 精度 的 地 理 信息 数 
据 。 随 着 3S 技 术 的 迅猛 发 展 ,传统 的 数据 获取 方法 


Rapid Mass Movement 


综 上 所 述 ,本 文 将 无 人 机 倾斜 摄影 技术 应 用 于 
新 疆 西 天 山 阿尔 先 沟 雪 朋 调查 人 研究 工作 中 ,利用 航 
拍 数 据 精 细 探 测 雪 崩 活 动 过 程 识别 雪崩 流 动 特征 
及 流 态 信息 ,生成 高 精度 DEM 数 据 ;再 结合 现场 实 
测 值 ,确定 RAMMS 模 型 输入 参数 ,进行 雪 朋 运动 过 
程 模拟 。 通 过 对 不 同类 型 雪崩 事件 的 模拟 和 重建 ， 
对 比分 析 传 统 地 面 调查 .无 人 机 遥感 解 译 结果 与 模 
拟 结果 的 差异 ,探讨 不 同类 型 ,不同 雪 层 释放 条 件 
下 雪 衣 活动 过 程 ,可 为 今后 雪崩 潜在 危险 预测 、 风 
险 规 避 及 灾害 应 急 处 置 提 供 强 有 力 的 数据 支撑 及 
科学 依据 。 


1 研究 区 概况 


天 山 山 脉 呈 准 东西 走向 ,地 形 构造 特殊 , 受 西 
风 带 气流 影响 ,在 西天 山 伊 犁 河谷 及 两 侧 山 区 易 形 
成 丰富 降水 区 ; 因 典 型 的 裙 皱 断 块 山 ,使 得 古 冰 斗 、 
雪 蚀 洼地 散 见 ,其 内 可 聚集 大 量 积 雪 , 因 而 成 为 雪 
表 的 多 发 地 2 2 。 而 阿尔 先 沟 正 处 于 西天 山 深 处 ， 
属 中 高 海拔 的 山地 地 貌 , 是 沿 纬 向 构造 山脉 发 育 的 
山 间 盆地 和 谷地 ,沟谷 两 侧 山 体 地 形 起 伏 较 大 , 山 
势 陡 峭 , 斜 坡 坡度 一 般 40°~60° ,地 面 高 程 2000 ~ 
3500 m; 河 谷 东 南 侧 林 间 少 有 雪崩 痕迹 ,而 西北 侧 
山地 裸露 ,植被 覆盖 率 低 ,属于 典型 的 雪崩 灾 害 孕 
灾 环 境 。2023 年 冬 春季 ,对 阿尔 先 沟 进行 雪 衣 调查 
发 现 ,该 段 在 14 km 的 主要 雪崩 区 内 ,发 育 雪 月 近 百 
处 , 且 发 生 频 繁 而 集中 。 阿 尔 先 沟 区 位 及 沟 内 实况 
如 图 1 所 示 。 

雪 情 很 大 程度 上 决定 着 雪 骨 触发 及 其 破坏 力 
的 大 小 31。 通过 建设 DZZ4 型 区 域 自动 气象 站 监测 
阿尔 先 沟 2022 一 2023 年 冬 春季 雪 情 变化 情况 ,对 于 
气温 、 风 向 等 变 率 较 大 的 观测 指标 ,采用 1 s 一 次 的 
采样 方法 ,风速 采用 0.25 s 一 次 的 采样 方法 ; 积 雪 深 
BE .固态 降水 等 采用 一 分 钟 一 次 的 采样 方法 。 气 象 


存在 采集 率 低 .生产 成 本 高 . 建 模 周期 长 .精度 低 等 
问题 ,不 能 满足 雪 骨 重建 需要 ” 。 而 利用 无 人 机 进 
行 数 据 采集 实景 建 模 等 技术 日 趋 成 熟 , 其 中 ,无 人 
机 倾斜 摄影 技术 可 通过 多 镜头 的 倾斜 相机 采用 定 
点 上 曝光 获取 地 面 多 角度 纹理 影像 ", 再 通过 处 理 生 
成 实景 三 维 模型 ,在 此 基础 上 生成 的 数字 高 程 模型 
(Digital Elevation Model, DEM) 及 数字 正 射 影像 
(Digital Orthophoto Map, DOM ) 精 度 高 .成 像 清 晰 ,是 
雪崩 运 动 过 程 模 拟 中 可 靠 的 数据 源 。 


站 具体 位 置 及 情况 如 图 1 所 示 。2022 一 2023 年 阿尔 
先 沟 冬 季 积 雪 深 厚 , 旦 持续 时 间 长 ,降雪 期 从 9 月 一 
次 年 4 月 结束 ,而 11 月 开始 平均 气温 低 于 0% ,地 表 
开始 积 雪 ,地 表 平 均 积 雪 深 度 最 深 可 达 91 om; 11 月 
一 次 年 4 月 ,气温 存在 明显 的 回升 情况 ,同时 河谷 内 
出 现 盛行 风 , 且 在 4 月 出 现 积 雪 期 内 平均 风速 的 峰 
值 点 ,达到 2.87 m.s ,其 最 大 风速 可 达 12 m-s” CS] 
2)。 由 于 风速 的 变化 使 风 吹 雪 形成 的 雪 层 在 密度 
和 硬度 上 都 呈现 很 大 的 变化 , 雪 层 中 的 应 力 在 局 部 
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(c) 阿尔 先 沟 冬季 实况 


(b) 阿尔 先 沟 区 域 位 置 
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(d) 阿尔 先 沟 气象 站 


图 1 阿尔 先 沟 地 理 位置 及 实况 图 


Fig. 1 Location ofAerxiangou and the actual map 
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(b) 风速 及 地 表 平 均 积 雪 深 度 变化 
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图 2 阿尔 先 沟 气象 站 2022-01- 一 2023-04 月 气象 监测 数据 


Fig.2 Meteorological monitoring data between 2022-01 and 2023-04 at Aerxiangou meteorological station 
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过 程 中 新 降雪 密度 的 增加 ( 即 高 密度 降雪 位 于 低 密 
度 雪 层 之 上 ) ,导致 雪 层 不 稳定 性 增 大 ,在 风力 和 降 
雪 等 外 力作 用 下 ,研究 区 雪崩 爆发 概率 增 大 ”。 


2 数据 处 理 与 方法 


21 空 -地 协同 调查 
利用 人 工地 面 调 查 方法 实地 踏勘 雪 朋 灾害 现 


场 , 获 取 堆 积 区 范围 .形态 及 积 雪 特征 等 实测 数 
据 。 联 合 以 无 人 机 航拍 为 核心 的 雪崩 无 人 机 空中 
调查 技术 获取 高 精度 航拍 数据 ,结合 航拍 数据 建立 
三 维 模型 ,对 调查 区 域内 的 雪崩 灾害 点 和 隐患 点 进 
行 系 统 性 全 面 分 析 , 判 断 雪崩 类 型 ,识别 雪崩 释放 
区 运动 区 及 堆积 区 , 查 明 灾 害 空间 分 布 状态 及 数 
量 。 在 此 基础 上 进行 雪崩 过 程 数 值 模 拟 , 对 比分 析 
堆积 体 与 野外 观测 范围 ,以 评价 模拟 结果 的 准确 
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性 ;通过 三 维 模型 解 译 .数值 模拟 结果 与 现场 调查 
数据 相互 验证 ,完善 灾害 点 基本 属性 信息 ,构建 雪 
骨 灾 害 数据 库 , 为 灾害 防治 提供 依据 。 具 体 雪 骨 灾 
害 空 -地 协同 调查 流程 如 图 3 所 示 。 


LT LEAT 
倾斜 摄影 航拍 数据 


利用 Agisoft Metashape "v 
Professional v 4 prs a r i 


多 视点 三 维 模型 
雪 骨 运动 特性 分 析 


雪崩 过 程 数 值 模拟 


图 3 空 -地 协同 调查 流程 图 


Fig.3 Air-to-ground cooperative investigation flow chart 


21.1 地 面 调查 在 传统 地 面 调查 过 程 中 对 调查 区 
域 按 照 雪崩 类 型 进行 分 区 判断 。 最 终 确定 重点 调 
查 区 域 一 般 调 查 区域 非 调 查 区 域 。 重 点 调查 区 
域 是 指 该 区 域 满 足 雪 崩 灾 害 诱发 条 件 旦 对 公路 、 重 
要 构造 物 等 产生 较 大 影响 。 非 调查 区 域 是 指 该 区 
域内 雪崩 灾害 不 会 对 公路 重要 构造 物 等 产生 影 
响 。 一 般 调查 区 域 则 是 指 雪崩 灾害 诱发 条 件 较 差 ， 
可 能 对 公路 、 重 要 构造 物 等 产生 较 小 影响 (图 4)。 


Y v 
一 般 调查 区 域 | [3E 81 Db 


一 般 调 查 区 域 


重点 调查 区 域 
无 人 机 雪崩 调查 


图 4 雪 骨 灾害 分 区 判断 流程 图 
Fig. 4 Avalanche disaster zoning judgment flowchart 


针对 重点 调查 区 域 进 行 现场 测量 获取 雪 骨 灾害 实 
测 数据 并 采用 无 人 机 雪 朋 调查 技术 进一步 解 译 雪 
RE o 

2.1.2 无 人 机 空中 调查 及 数据 获取 根据 地 面 调查 
分 区 结果 ,在 重点 调查 区 域 开展 无 人 机 雪崩 调查 。 
前 期 准备 工作 包括 收集 气象 数据 .定点 勘察 和 划 定 
飞行 区 域 ;在 此 基础 上 规划 航线 ,设置 相关 飞行 参 
数 。 为 保证 测 区 完整 性 设置 了 10%~30% 的 缓冲 区 ， 
同时 为 确保 无 人 机 航拍 能 够 尽 可 能 多 的 覆盖 测 区 
范围 ,根据 测 区 及 缓冲 区 大 小 设置 45° 的 最 佳 倾 斜 摄 
影 角 度 。 通 过 提高 相片 重 登 度 及 拍摄 间隔 方式 以 保 
证 测 区 范围 内 各 处 至 少 有 三 幅 不 同位 置 的 重 秋 影像 
保证 了 模型 细节 的 清晰 以 及 数据 的 精确 性 ”你 。 结 
合 无 人 机 倾斜 摄影 三 维 重 建 技术 从 一 个 垂直 ,多 个 
倾斜 视角 出 发 ,获取 不 同 地 表 地 物 顶 部 信息 ,以 及 
立 面 详细 的 轮廓 和 纹理 信息 0907" ,识别 雪 骨 释放 
区 运动 区 及 堆积 区 , 查 明 雪崩 空间 分 布 状态 。 在 
三 维 重 建 过 程 中 可 对 航拍 数据 进行 预 处 理 以 解决 
图 像 失 真 .变形 等 问题 ;处 理 方式 包括 利用 相机 畸 
变 参 数 对 图 像 畸 形 矫 正 以 及 考虑 到 航拍 过 程 中 由 
于 天 气 原因 造成 的 曝光 或 色彩 不 均 问 题 进行 的 羽 
化 及 重 曝 光 人 处理 淖 。 经 过 预 处 理 后 的 航拍 数据 导 
入 多 视点 三 维 建 模 软 件 获取 密集 点 云 数据 ,生成 三 
维 模 型 ,导出 DEM, DOM, SERA AY Be RE AT Jo 
ABL FEAT , Dr AS [a] HE. TEA 
DLE AI ep CS LEE AL 5 


in 


前 期 准备 


|>| 航线 规划 Hirr] 


«——] 


获取 密集 点 云 数 据 航拍 数据 处 理 


DOM | 


| = 三维 模型 | | DEM | 
[ 


遥感 解 译 


图 5 无 人 机 雪崩 调查 技术 流程 图 
Fig. 5 Flow chart of UAV avalanche survey technique 


2. RAMMS 雪 骨 过 程 模 拟 原理 及 输入 数据 获取 
2.2.1 RAMMS 模型 模拟 原理 分 析 采用 RAMMS 模 
型 ,模拟 复杂 地 形 中 雪崩 的 流动 过 程 ,其 本 质 是 在 
VS 模型 基础 上 考虑 雪崩 运动 过 程 中 流动 高 度 和 速度 
变化 的 非 定 值 和 非 均匀 运动 特征 ,使 用 积 雪 深 度 平 
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均 方 程 同 时 耦合 RKE 模 型 来 考虑 质量 和 动能 守恒 ， 
解释 雪 粒 之 间 的 随机 运动 和 非 弹性 互相 作用 。 
模型 基于 笛 卡 尔 坐 标 系 来 表征 平均 速度 U0, 假 
设 和 和 Y 为 固定 于 笛 卡 尔 坐标 系 中 的 水 平 坐标 ,Z 
(x,y) 表 示 用 XX 和 了 参数 化 的 山脉 轮廓 , 自 变 量 X、Y 
表示 地 形 弧 长 ,坐标 XY.Z 形 成 了 曲面 诱导 坐标 
系 。 随 着 时 间 1 的 变化 雪 朋 流动 高 度 H(x,y,1)(m) 
和 平均 速度 U(x,y,t) (ms ) 也 随 之 变化 ,平均 速度 
U(x,y,t) 关 系 式 如 下 : 
Ux,y,d)= [U, (.y. Dy] (1) 
n, = HU., Uy, (2) 


式 中 : 忆 为 X 轴 的 流速 值 ; 忆 为 了 轴 的 流速 值 ;7 为 平 
均 流 速 矩 阵 的 转 置 妹 表 示 雪 衣 运 动 的 方向 。 
由 于 雪崩 呈现 浅 层 流 体 几 何 形态 ,用 质量 平衡 
方程 来 求解 雪崩 流 动 高 度 万 ,如 下 式 : 
0,+90,(HU,)+0,(HU,)=Q(y,d) G) 
VS 模型 中 将 摩擦 系数 考虑 为 两 部 分 ,与 法 回应 
力 N 成 正比 的 干 库伦 摩擦 系数 i 和 满 流 摩擦 系数 &。 
摩擦 阻力 S(Pa) 公 式 如 下 : 
S=uN+ wp N=phe cos (4) 
SUH: p 表示 密度 ;g 为 重力 加 速度 ;4 为 斜 角 ; h Nit 
BPE RE su NCE um (u, su, ;摩擦 系数 对 碎 届 流 运 
动 演变 过 程 有 直接 影响 。 当 流动 接近 于 停止 时 , 
起 主导 作用 ; 当 流 速 快速 运动 时 ,é 主 导 作用 。 
考虑 到 雪崩 运动 体 是 由 雪 团 组 组 合 , 运 动 过 程 
中 可 能 携带 粒状 碎 导 等 物质 ,运动 过 程 中 平行 或 者 
非 平 行 斜坡 方向 时 会 出 现 速 度 上 的 波动 ,采用 RKE 
模型 实时 修整 ,可 更 好 的 表现 雪崩 实 时 运动 特征 的 
变化 。 如 公式 (5) 和 公式 (6) 所 示 : 


uR)-2u ex(- E (5) 


TERES (6) 


式 中 :Ro 为 平均 随机 动能 密度 函数 的 摩擦 指数 增长 
率 , 为 常数 ;RR 为 深度 平均 随机 动能 。 修 正 后 的 干 库 


伦 摩擦 系数 4 和 漠 流 摩擦 系数 E 取 决 于 R。 

2.2.2 雪 层 断裂 深度 的 确定 ”断裂 深度 是 指 雪 骨 释 
放 区 内 ,一 定 厚度 的 雪 层 在 重力 温度 等 因素 作用 
下 , 积 雪 薄弱 层 失 效 出 现 垂直 于 坡 面 的 裂缝 后 , 沿 
着 裂缝 向 下 滑动 的 平均 积 雪 厚 度 。 根 据 山 坡 积 雪 
体 的 平衡 条 件 , 山坡 上 积 雪 所 受 重力 可 分 解 为 平行 
于 山坡 的 应 力 和 垂直 于 山坡 的 正 压力 ,而 下 雪 层 在 
重力 的 作用 下 有 沿 坡 面 向 下 滑动 的 趋势 , 当 雪 层 的 
下 滑 力 大 于 阻力 时 引发 雪崩 。 由 此 得 出 山坡 坡 面 
积 雪 临界 厚度 ,如 公式 (7) 所 示 。 当 坡 面 积 雪 厚度 
到 达 临 界 厚 度 值 时 , 雪 层 的 下 滑 力 与 阻力 处 于 平衡 
状态 ,这 时 如 果 积 雪 厚 度 的 增 大 或 者 雪 层 应 力 分 布 
ARIS] ,都 会 导致 雪 层 沉降 ,引起 滑 塌 。 因 此 ,本 研究 
引用 临界 厚度 来 确定 雪 层 断裂 深度 。 


- c 
De a eaei () 


SU: h, Ay LU AR JE E ; c Ay Fy 53 3C TB TRE BE BR 
J; p HE 5j Hi IC BS LEE HE f (9) ;7 为 积 雪 的 密度 
(genm”);a 为 山坡 坡度 (?)。 

由 于 公式 (7) 中 的 变量 均 与 积 雪 、 坡 面 特性 相 
关 ; 结 合 研 究 区 现场 实际 情况 与 胡 汝 双 等 ”3 对 天 
山 山 区 季节 性 积 雪 物理 性 质 研 究 过 程 中 得 出 的 不 
同类 型 积 雪 物理 力学 强度 分 析 结 果 ,确定 公式 (7) 
中 数值 取 值 ( 表 1) ,其 中 山坡 坡度 根据 地 形 分 析 
得 出 。 
2.2.3 摩擦 系数 的 确定 RAMMS 对 摩擦 系数 (4) 
有 两 种 设置 情况 ,常量 和 变量 计算 模式 。 常 量 计算 
模式 考虑 雪 朋 区 域 下 垫 面 处 于 相同 性 质 的 地 表 状 
DU ,整个 雪 朋 运动 区 域 (释放 区 、 运 动 区 堆积 区 ) 用 
固定 的 wc 和 # 值 ,摩擦 系 数 (w 6) 的 取 值 取决 于 重 现 期 
以 及 雪 骨 体积。 而 变量 计算 模式 ,yw 和 é 值 不 固定 ， 
考虑 不 同 海拔 区 域 地 表 状 况 不 同 ,使 用 不 同 的 wkc 和 & 
值 。 根 据 研究 区 内 雪 前 特征 并 结合 模型 参数 建议 
值 进 行 匹配 以 适用 于 研究 区 域 。 
2.2.4 数字 高 程 模型 DEM 作为 RAMMS 数值 模拟 
中 最 基本 的 输入 数据 ,能 够 准确 反映 雪崩 模 拟 区 域 
的 地 形 特征 ,而 其 分 辩 率 的 高 低 对 于 模拟 效果 存在 
BETZ MOO” | Christen?) 分别 采用 2.5 m.25 mm 计算 


表 1 积 雪 物理 性 质 


Tab.1 Physical properties of snow 


雪 型 ARIE sem?) 内 摩擦 系数 


密度 /(g.cm”) 抗 断 强度 /(g* cm?) 


新 雪 5 0.22 


0.15( 实 测 ) 2.8 
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网 格 分 状 率 进行 模拟 ,结果 显示 不 同 计算 网 格 分 辨 率 
下 雪崩 流 中 心 堆积 存在 分 布 差异 。 同 时 Buhler 等 ” 提 
出 25 m 的 空间 分 辨 率 足够 准确 的 模拟 大 规模 雪崩 ， 
较 2 m DEM 的 模拟 结果 具有 和 较 高 的 吻合 度 。 本 研 
究 在 考虑 区 域 面积 .计算 时 间 与 雪崩 规模 的 基础 
上 ,为 了 保证 模拟 结果 的 准确 性 将 采用 无 人 机 倾斜 
摄影 航拍 数据 获取 的 2.0m B DEM 作为 输入 数据 ， 
进行 雪 衣 过 程 模拟 。 
2.3 三 维 计算 模型 建立 

依据 无 人 机 倾斜 摄影 技术 获取 的 高 精度 DEM， 
通过 ArcMap 完成 地 形 表面 分 析 自 动 识 别 坡度 适 
合 区 、 划 定 雪崩 释放 区 域 以 及 森林 、 基 础 设施 等 位 
置 ” ,具体 程序 如 下 :(1) 创建 地 图 ,导入 DEM 创 建 
2D 环境 模型 ,(2) 执行 “创建 Shapefile” ,要素 类 型 为 
“Thi” ,坐标 系统 一 为 WGS_1984 ,根据 自动 识别 坡度 
结果 划 定 雪 骨 释放 区 域 ,并 释 加 DOM; 由 于 RAMMS 
对 于 摩擦 系数 的 设 定 分 为 常量 和 变量 计算 模式 ,如 
果 使 用 恒定 摩擦 值 进 行 计算 则 不 考虑 地 形 起 伏 和 
森林 区 域 ,而 使 用 变量 摩擦 值 则 需要 根据 地 形 数据 
对 摩擦 值 进 行 分 类 设置 。 因 此 , 需 对 森林 、 基 础 设 
施 等 区 域 进 行 划 分 以 便于 摩擦 值 分 类 设置 .(3) 将 
处 理 后 的 DEM、DOM 、Shapefile 文 件 导 入 雪崩 动 力 
学 数值 模拟 软件 一 RAMMS 中 ,再 根据 前 期 现场 调 
查 和 实测 数据 确定 模拟 软件 输入 参数 值 。 最 后 输 
出 雪 骨 三 维 计算 模型 。 


3 阿尔 先 沟 雪 毅 调 查 及 过 程 模 拟 


3.1 雪 朋 信息 提取 

2023 年 2 月 11 一 18 日 ,调查 发 现 A、B、C 三 处 雪 
前 ,根据 现场 信息 提取 , 雪 前 发 生 于 东北 -西南 向 
Y029 公 路 的 北 侧 陡 峭 山 体 ,岩石 裸露 e TA RS 
动 痕迹 明显 。 雪 骨 从 坡 面 陡峭 处 释放 , 冲 出 公 
制 动 区 位 于 公路 南 侧 硬 路 肩 外 。 对 A、B、C 三 处 雪 
崩 堆 积 区 进行 分 析 发 现 ( 图 6) ,A 处 雪崩 堆 积 体 剖 面 
结构 中 雪 衣 碎片 粒 径 较 小 且 以 粉 状 为 主 ,剖面 积 雪 
结构 较 均 匀 。 山 坡 坡 面 有 风 吹 雪 痕迹 。 经 现场 测 
量 ,A 处 雪 骨 堆积 区 抛 出 距离 约 占 公路 2 m, 堆 积 体 
约 为 42.83 mm ,高 度 约 为 1.7 m。 从 现场 情况 看 雪崩 


值 ,降雪 及 风力 等 看 合作 用 下 产生 的 新 雪 雪 月 。 

分 析 B 处 雪 骨 堆积 区 后 发 现 ,在 均匀 积 雪 结构 
堆积 体 的 上 方 有 少量 颗粒 较 小 的 雪崩 碎片 ;由 于 温 
度 上 升 , 雪 层 融化 坡 面 呈现 出 鱼 磷 状 雪 层 (图 6c)。 
依据 王建 等 提出 的 中 国 积 雪 特 性 现场 测量 方法 
进行 积 雪 特性 测量 中 发 现 , 坡 面 表层 颗粒 雪 密 度 达 
到 434 gcmi ,表层 温度 1.9 °C. MSE RK PS 
层 温度 呈现 出 逐步 增 大 的 趋势 ,地 表面 积 雪 层 100 
cm 和 0 cm 处 温度 分 别 为 -0.1% .0s%。 雪 层 密度 随 
着 温度 的 增 大 不 断 增 大 ,100 cm .0 cm 处 雪 层 密度 
分 别 为 263 g*cm?,385 gem? C 2, Al 6c) ,. TR ME 
积 体积 雪 密 度 的 变化 以 及 地 面 调查 情况 分 析 ,B 处 
堆积 区 是 在 第 一 次 新 雪 雪 前 后 形成 的 堆积 形态 基 
础 上 ,受到 温度 升 高 坡 面 上 方 积 雪 融化 后 ,由 近期 
再 次 出 现 的 小 规模 湿 雪 雪崩 形 成 的 堆积 体 。 

与 A.B 处 雪崩 不 同 ,C 处 属于 连续 性 雪崩 点 , 堆 
积 区 长 度 达 150 mm 左右 ,形态 呈 锥 型 (网 6d)。 堆 积 
区 截面 最 大 雪 深 1.8 m ,平均 雪 深 1.7m, 抛 出 距离 最 
远 约 占 公 路 8 m 宽 。 根 据 现 场 情况 分 析 , 山 体 下 方 
存在 断面 ,断面 上 设置 防护 网 ;山体 上 方 存在 类 似 
沟 体 的 结构 。 就 积 雪 密度 而 言 C 处 较 A、B 处 小 , 雪 
体 松散 ,未 出 现 冰 体 ,堆积 体 以 粉 状 雪 为 主 , 结 构 较 
均匀 , 属 新 雪 雪 前 。 

3.2 模型 构建 及 灾害 解 译 
3.2.1 影像 获取 本 次 调查 采用 大 吐 经 纬 M300 
RTK 无 人 机 ,搭载 禅 思 H20 摄像 头 ,测绘 区 域 面 积 
8.9256x10* m? 航 线 时 长 约 为 32 min ,飞行 参数 如 表 3 
所 示 。 本 次 飞行 任务 共 5 条 航线 (图 7) ,获得 58 张 
倾斜 摄影 图 像 。 由 于 高 处 、 低 处 分 辩 率 不 同 ,其 影 
像 地 面 分 辨 率 为 0.05~0.10 mo 
3.2.2 数据 处 理 及 灾害 解 译 采用 Agisoft Metashape 
Professional 对 航拍 数据 进行 处 理 , 生 成 三 维 实景 模 
型 ,处理 步 又 包括 对 齐 照 片 .建立 密集 点 云 . 生 成 网 
格 .生成 纹理 ;根据 密集 点 云 数 据 生 成 正 射影 图 像 ， 
在 正 射 影 图 像 基 础 上 进行 三 维 实景 模型 建立 (图 8)。 
结合 地 面 调 查 情况 以 及 地 形 判 断 ,A、B 处 属于 
坡 面 型 雪崩 ,该 坡 面 上 方 存在 两 个 天 然 的 积 雪 平 
台 , 汇 雪 面 积 在 0.58239x10~1.3108x10* m^, HK Tj 


释放 区 可 能 位 于 位 置 较 低 且 面积 较 小 的 局 部 坡 面 ， 
雪崩 流 在 坡度 小 于 25° 的 位 置 停止 释放 (图 6b)。 初 
步调 查 判 断 该 区 域 是 由 于 坡 面积 雪 厚 度 达 到 临界 


呈现 出 上 陡 下 缓 的 形状 ,上 半 部 分 陡坡 处 坡度 达 
50" 以 上 , 坡 面积 雪 量 少 ,地 表 裸露 , 坡 面积 雪 主要 集 
中 在 下 半 部 分 缓坡 处 ,坡度 在 30*~45*; 积 雪 平台 至 
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(a) 雪崩 点 区 位 图 


84°40'30"E 84°41'0"E 84°41'30"E 
T T T 


43°23'0"N 


43°22'30"N 
43°22'30"N 


84°40'30"E 84°41'0"E 84°41'30"E 
(b) A 处 雪崩 堆 积 区 
| 


f 


af a 
(c) B 处 雪崩 堆积 区 


截面 最 大 雪 深 : 
1.8m 


图 6 雪崩 点 位 置 及 堆积 区 实况 图 


Fig.6 Actual map of avalanche point location and accumulation area 
坡 底 最 短 距离 在 510 m ,最 长 距离 为 594 m( 图 8)。 较 广 ,堆积 区 积 雪 量 较 大 ,但 整体 冲击 力 流动 速度 
造成 坡 面 型 雪 骨 发 生 的 主要 原因 :(1) 缓坡 坡 BE (2) 上 方 积 雪 平台 雪 层 失 稳 , 积 雪 释 放量 较 
面积 雪 达 到 临界 厚度 出 现 雪 层 失 稳 , 雪 骨 影 响 范 围 。 大 ,两 个 积 雪 平台 共用 同一 个 运动 区 , 易 造 成 两 个 


HH 
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表 2 堆积 区 积 雪 特性 


Tab.2 Snow cover characteristics in accumulation area 


雪 层 厚度 /em 
100 80 60 40 20 10 0 
雪 层 温度 /*C -0.1 -0.6 -0.4 -0.1 -0.1 0 0 
SEES BEI (g* em?) 263 283 303 314 337 360 385 


表 3 无 人 机 倾斜 摄影 飞行 参数 
Tab.3 UAV tilt photography flight parameters 


项 目 参数 项 目 参数 
模式 倾斜 摄影 | 拍照 模式 等 距 间 隔 拍照 
航线 /条 5 PIDE: E 75 
倾斜 角度 /(°) 45. [MREZE K% 80 
飞行 高 度 /m 1000 | 旁 向 重 赤 率 ( 倾 斜 )/% 70 
边 距 /m 20 “| 航线 重生 率 (倾斜 /9% 70 


图 7 研究 区 域 航线 图 
Fig.7 Study the regional route chart 


ELE S 
BAEHERRIX 710m 


TP, z 
C 处 堆积 区 ah 


积 雪 区 同时 释放 ,雪崩 流 具 有 交汇 点 。 考 虑 积 雪 夹 
带 作用 , 易 引 起 缓坡 坡 面 积 雪 的 失 稳 。 和 危险 性 较 
大 ,影响 范围 较 广 。 

结合 地 面 调查 情况 以 及 地 形 判 断 ,C 处 堆积 区 
属于 沟 模 一 坡 面 复合 型 雪 朋 ,上 部 释放 区 为 沟 权 发 
育 , 上 宽 下 窜 , 最 宽 处 大 于 300 m, 最 窄 处 不 足 50 m。 
汇 雪 面 多 , 积 雪 量 大 ,整体 汇 雪 面积 可 达 38124.0 m, 
释放 区 至 坡 低 直 线 距 离 近 700 m ,落差 较 大 。 经 测 
量 汇 雪 区 山体 坡度 在 26°~35° ,具有 较 好 的 储 雪 条 
件 。 雪 崩 的 发 生 主 要 由 于 沟 槽 型 山体 左右 两 侧 坡 
面积 雪 失 稳 汇 集 于 沟 槽 内 ,上 宽 的 沟 体 使 得 汇 雪 量 
大 ,下 案 使 得 雪崩 冲 出 后 到 达 坡 面 旦 发 散 状 向 左右 
两 侧 运 动 (图 8, 图 9)。 

根据 航拍 数据 解 译 出 该 区 域 山体 存在 雪 榴 和 
雪 层 裂缝 。 利 用 无 人 机 Real Time Kinematic (RTK ) 
空中 测量 技术 发 现 坡 面 型 雪崩 中 以 两 条 雪 构 规模 
最 大 ,长 度 在 80~150 mW LUA A AI ak, SEE 
走向 。 雪 柳 旁 伴随 雪 层 裂 颖 发育, 裂缝 较 长 ,两 条 
裂缝 长 约 155.41 m、144 m, 裂 颖 周围 有 较 明显 的 滑 
动 痕迹 。 下 半 部 分 坡 面 处 出 现 竖 直 向 下 的 裂缝 , 裂 


T 


三 维 实景 模型 
3D real life model 
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图 9 雪 骨 解 译 结果 


Fig.9 Disaster interpretation results 


颖 较 长 ,达到 254.28 mm ,延伸 至 坡 底 (图 9c)。 对 沟 酸 
一 坡 面 复合 型 雪崩 的 灾害 解 译 过 程 中 发 现 该 雪 朋 
汇 雪 区 出 现 多 个 雪 柳 ,规模 最 大 的 长 度 为 191.88 m, 
呈 东 西 一 南北 走向 , 雪 层 裂 颖 与 雪 榴 平行 发 育 ,长 
度 达 184.38 m。 释 放 区 域 可 见 多 条 裂 颖 以 及 雪 衣 请 
动 痕迹 (图 9)。 

结合 解 译 结果 与 地 面 调查 综合 分 析 , 目前 研究 
区 内 仍 处 于 灾害 发 育 阶段 ,如 果 持 续 降 雪 后 温度 持 
续 上 升 , 雪 榴 裂缝 具有 加 剧变 形 的 趋势 ,将 会 出 现 
湿 雪 型 雪 板 雪崩 。 雪 层 中 含水 量 较 高 ,导致 衣 出 距 
离 大 ,从 而 对 环境 和 基础 设施 更 具 威 胁 。 其 次 ,还 
易 在 崩落 过 程 中 庄 挟 土壤 ,造成 土壤 流失 3。 
3.3 雪 骨 过 程 模拟 及 分 析 

将 无 人 机 航 摄 数据 处 理 生 成 的 DEM 通过 Arc- 
Map 完成 前 期 工作 后 ,结合 地 面 调查 以 及 无 人 机 航 
拍 数据 ,确定 模型 输入 参数 ,再 导 人 RAMMS 对 区 域 
内 坡 面 型 , 沟 槽 一 坡 面 复合 型 雪崩 进行 模型 计算 ， 
分 析 灾 害 活动 特征 。 本 文 考虑 最 大 危险 性 , 即 雪 骨 
量 最 大 的 情况 (山坡 积 雪 全 层 断 裂 ) ,认为 临界 厚度 
等 于 断裂 深度 。 


坡 面 型 h = C 
Tm Te Tr(sina - cos aX tgo) 
5 
= = 0.6 
0.15(sin 45°- cos 45°X0.22) ~" 
(8) 
HK: h,= C 
m : T(sina -cosa X ted) 
i 5 
0.15(sin 35? — cos 35? x 0.22) 
£0.84 m 
(9) 


考虑 到 本 文人 研究 的 坡 面 型 雪 骨 属于 非 沟 洪 地 
形 , 且 雪 骨 运动 区 域 (释放 区 运动 区 .堆积 区 ) 的 下 


垫 面 (剥蚀 构造 山区 ) 一 致 ,摩擦 系数 采用 常量 计算 
模式 ,而 沟 档 一 坡 面 型 雪崩 流 经 过 区 域 为 沟渠 \ 沟 
谷山 谷 平坦 地 区 ,摩擦 系数 则 采用 变量 计算 模 
式 。 因 此 本 次 雪崩 模拟 所 需 的 摩擦 系数 需要 结合 
地 形 特 性 和 模型 参考 值 进行 分 类 选择 。 具 体 摩擦 
系数 参考 值 如 表 4 所 示 。 


表 4 摩擦 系数 4 和 c 在 不 同 地 形 \ 不 同 海拔 高 度 取 值 
Tab.4 The values of friction coefficients u and ¢ for 


different terrain and altitudes 


地 形 条 件 海拔 /m m É 

非 沟渠 型 地 形 >3000 0.29 1500 
3000~2700 0.305 1400 

<2700 0.33 1250 

沟渠 型 地 形 > 3000 0.34 1250 
3000~2700 0.355 1180 

«2700 0.37 1050 

沟谷 地 形 >3000 0.45 1050 
3000~2700 0.46 1100 

<2700 0.47 900 

平坦 地 形 >3000 0.275 1750 
3000~2700 0.285 1600 

<2700 0.295 1500 


3.3.1 坡 面 型 雪崩 模拟 结果 及 分 析 根据 模拟 结 
显示 ,以 坡 面 上 方 积 雪 平台 为 释放 区 , 雪 骨 释放 体 
积 可 达 8.2669x10 m。 雪 月 运动 方向 主要 沿 坡 面 运 
动 ,雪崩 流 较 集 中 ,支流 少 。 堆 积 区 分 布 集中 ,影响 
范围 较 集中 。 两 个 积 雪 平台 共用 同一 个 运动 区 , 相 
互 影 响 下 连续 性 释放 , 雪 骨 流 具 有 交汇 点 ,与 前 文 
三 维 实景 模型 解 译 结果 相同 。 但 产生 的 堆积 区 分 
布 集中 , 且 连 续 ,与 地 面 调 查 过 程 中 发 现 的 分 散 的 
堆积 区 布局 不 同 。 进 而 判断 近期 坡 面 型 雪崩 并 非 
职 雪 平台 的 释放 。 
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以 积 雪 平 台 为 潜在 释放 区 考虑 雪 骨 量 最 大 时 
( 雪 层 全 层 断 裂 ) ,雪崩 流 运动 时 长 仅 为 128 s, 并 且 
40s 内 雪 朋 流 已 到 达 公 路 , 制 动 区 位 于 公路 南 侧 硬 
路 肩 外。 依据 地 形 分 布 堆积 体 呈 现 出 中 间 高 ,四周 
低 的 椎 状 体 。 由 图 10 可 知 ,堆积 区 约 在 120s 内 流 
动 高 度 达 到 最 大 , 约 为 3.55 m; 而 最 大 流动 速度 在 
18.34 m-s”, 冲击 力 最 大 在 32.67 kPa; 积 雪 堆 积 量 约 
为 1.8525x10' ,影响 面积 在 3369.7 m^ 左右。 说 明 
该 情况 下 雪 骨 影响 范围 广 ,堆积 体面 积 较 大 ,冲击 
力 及 运动 速度 足以 对 环境 和 基础 设施 等 承载 体 产 
生 较 大 威胁 。 
3.32 沟 楼 一 坡 面 复合 型 雪 前 模拟 结果 及 分 析 E 
虑 雪崩 量 最 大 的 情况 ( 雪 层 全 层 断 裂 ) 进 行 模拟 , 沟 
柳 一 坡 面 复合 型 雪 朋 体 从 海拔 2988 m 人 处 启动 ,释放 
体积 可 达 7.6071x10' m: ,从 坡 面 骨 塌 而 下 , 沿 着 上 宽 
下 罕 且 深 的 沟 体 流出 后 顺 着 坡 面 呈 扇 形状 流出 ,大 
约 在 75s 后 开始 堆积 在 海拔 约 2229 mm 平坦 区 ,256 s 
后 停止 堆积 ,形成 面积 为 1369.7 my, 体积 约 11717 m° 


(a) 雪 骨 流动 高 度 
84°40'30"E 
d XX 


E T 


84°41'0"E 


BEES 
一 堆积 区 
=y029 


2 雪崩 堆积 体 厚度 /m 
[m 
0 


84?41'0"E 


43?22'30"N 


84?40'30"E 


(c) 雪崩 冲 击 力 


43°22'30"N 


的 堆积 体 。 堆 积 体 最 大 堆积 高 度 位 于 河谷 内 ,可 达 
3.8 m, 平 均 堆 积 深 度 约 为 2.5 m, 流 动 速度 最 大 值 为 
13.6 m: s ,冲击 力 最 大 值 为 30.72 kPa。 由 于 沟 体 深 
MÆ ,运动 区 域 小 ,大 部 分 积 雪 遗留 在 沟谷 中 ,最 大 
堆积 深度 达到 18.91 m, 是 堆积 区 最 大 深度 的 两 倍 
多 。 并 且 该 雪崩 运动 过 程 中 各 运动 特征 高 值 均 位 
于 沟谷 中 ,最 大 速度 32.06 ms', 最 大 冲击 力 154.19 
kPa( 图 11a~ 图 11c)。 

沿 着 Y029 公 路 获取 堆积 体 横 截 面 , 根 据 横 截面 
堆积 高 度 反 应 出 堆积 体形 态 , 受 坡 面 形 态 影 响 形成 
中 间 高 .四 周 低 的 连续 性 椎 状 堆积 体 , 与 前 文 地 面 
调查 堆积 形态 和 航 摄 成 果 解 译 结果 一 致 (图 11j)。 
但 雪 层 全 层 断 裂 下 堆积 区 平均 堆积 高 度 是 现场 测 
量 的 堆积 区 平均 截面 深度 1.71 m 的 近 1.5 倍 。 与 实 
际 调查 堆积 情况 不 符 。 其 原因 与 释放 量 有 关 , 雪 骨 
模拟 过 程 中 堆积 区 雪 量 受到 摩擦 系数 和 释放 量 的 
控制 ,在 下 垫 面 处 于 相同 性 质 的 地 表 状 况 下 ,摩擦 
系数 相同 ,而 不 同 的 释放 量 对 雪 骨 运动 的 冲 出 范围 


(b) 堆积 区 积 雪 流动 高 度 变化 
40 —  VESUN HERE FE TR 


流动 速度 /ms 
25.56 


图 10 坡 面 型 雪崩 运动 特征 分 析 


Fig. 10 Analysis of slope avalanche movement characteristics 


11 期 张 天 意 等 :基于 空 -地 协同 调查 的 西天 山 阿 尔 先 沟 雪 骨 过 程 数值 模拟 1739 


(a) 雪 层 全 层 断 裂 下 雪崩 流动 高 度 (b) 雪 层 全 层 断 裂 下 雪崩 冲 击 力 (c) 雪 层 全 层 断 裂 下 雪崩 流 动 速度 
流动 高 度 /m 冲击 力 /kPa 流动 速度 /ms 


18.91 154.19 32.06 
i 9.46 so i 71.09 ü 16.03 
0 j 0 0 


(d) 雪 层 断裂 厚度 减少 20% (e) 雪 层 断裂 厚度 减少 20% O 雪 层 断裂 厚度 减少 20% 
下 雪山 流动 高 度 下 雪山 冲击 力 下 雪山 流动 速度 
流动 高 度 /m 冲击 力 /kPa 


» 16.62 3 139.55 
T. "A x "T 
0 0 


流动 速度 /(ms ) 


d 30.05 
ss ü 15.25 
T 0 


(g) 雪 层 断裂 厚度 减少 40% (b) 雪 层 断裂 厚度 减少 40% 
Hi PSA 


雪 层 断裂 厚度 减少 40% 
下 雪 衣 流动 高 度 mos 
流动 高 度 /m 


(OJE 
下 雪崩 流动 速度 


流动 速度 /ms 


18.41 
i 14.21 
0 
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: de 
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图 11 DEKRA NE as SET 


Fig. 11 Characterisation of trench-slope composite avalanche kinematic signature map 


1740 


ju 


4035 


及 堆积 量 产生 较 大 影响 。 

为 了 使 得 堆积 区 模拟 结果 更 符合 地 面 调查 情 
DL ,利用 单 因 素 分 析 方 法 进行 不 同 断 裂 深 度 下 雪崩 
数值 模拟 ,将 断裂 深度 分 别 减少 20% 40% ,所 得 的 
断裂 深度 分 别 为 0.67 m、0.34 m。 从 模拟 结果 表 5 可 
以 看 出 , 当 雪 层 断 裂 深度 减少 20%, 雪 月 流 积 区 影响 
范围 为 1299.5 m? ,平均 堆积 深度 为 2.27 m, 冲 出 量 为 
4578.22 m? ,堆积 区 最 大 流动 速度 为 10.24 m- s~, ie 
大 冲击 力 为 21.42 kPa; 当 雪 层 断 裂 深度 减少 40%， 
雪崩 堆积 区 影响 范围 在 1178.5 mv 左右 ,平均 堆积 深 
度 为 1.64 m, 冲 出 量 为 3107.76 mi’, 最 大 流动 速度 为 
6.58 ms ,最 大 堆积 区 冲击 力 为 17.97 kPa。 由 于 三 


破坏 性 损伤 ,严重 情况 下 可 威胁 生命 财产 安全 。 


4 结 论 

本 文 基于 空 -地 协同 调查 方案 对 雪崩 过 程 进行 
数值 模拟 并 对 运动 特征 进行 分 析 , 得 出 以 下 结论 ; 

(1) 以 无 人 机 倾斜 摄影 为 核心 的 雪 衣 灾 害 空 
中 、 地 面 协同 调查 ,将 传统 的 地 面 调查 方法 结合 无 
人 机 遥感 解 译 和 数值 模拟 结果 相互 验证 ,提高 了 对 
于 雪崩 灾害 发 育 状况 评价 的 准确 性 ,弥补 了 过 多 于 
依赖 经 验 值 判 断 的 主观 性 。 

(2) 2 月 中 旬 , 西 天山 阿尔 先 沟 为 持续 降雪 条 件 


次 模拟 在 雪 骨 释放 量 上 存在 很 大 的 差异 , 相 比 于 雪 
层 全 层 断 裂 , 断 裂 深度 减少 40% 雪 崩 流 到 达 堆 积 区 
时 间 较 晚 ,但 持续 时 间 最 长 ,接近 于 300 s ,而 断裂 深 
度 减 少 20% 雪 前 流 在 80 s 后 流 到 堆积 区 ,持续 时 间 
在 275 ss 左右, 持续 时 间 均 大 于 雪 层 全 层 断 裂 ;堆积 
区 积 雪 厚度 最 大 值 则 因 持 续 时 间 的 增长 形成 较 晚 
(图 11k, 表 5), 其 原因 主要 与 雪崩 流 自重 相关 。 


表 5 不 同 断裂 深度 下 模拟 结果 统计 
Tab.5 Statistics of simulation results at 
different fracture depths 


不 同 断 裂 深 度 /m 

0.84 0.67 0.34 
运动 时 长 /s 256 275 298 
堆积 量 /m 11717 457822 3107.76 
堆积 面积 /mY 1369.7 1299.5 1178.5 
堆积 区 最 大 流动 高 度 /m 9.22 7.03 4.23 
最 大 堆积 高 度 /m 3.8 3.05 2.47 
平均 堆积 深度 /m 2:5 220 1.64 
堆积 区 最 大 流动 速度 /ms ) 13.6 10.24 6.58 
堆积 区 最 大 冲击 力作 Pa 30.72 21.42 17.97 


由 以 上 结果 可 知 , 沟 枝 一 坡 面 复合 型 雪 骨 由 于 
地 形 因 素 影响 ,堆积 区 分 布 状况 主要 受 下 部 坡 面 控 
thi] ,在 对 运动 特征 进行 分 析 中 发 现 ,断裂 深度 减 小 ， 
导致 释放 量 减少 ,堆积 量 . 最 大 堆积 高 度 .最 大 堆积 
区 流动 速度 .最 大 堆积 区 冲击 力 随 之 减 小 , 而 堆积 
区 面积 整体 变化 幅度 较 小 。 观 测 值 与 模拟 结果 对 
比 发 现 ,断裂 深度 减少 40% ,平均 堆积 深度 在 1.64 m 
左右 ,与 现场 堆积 体 测 量 结果 相近 ,堆积 区 最 大 流 
动 速度 达到 了 6.58 ms ,堆积 区 最 大 冲击 力 在 
17.97 kPa 左右 ,最 大 运动 特征 值 足 以 对 承载 体 产生 


下 坡 面积 雪 厚 度 趋 近 于 临界 厚度 值 , 雪 层 失 稳 出 现 
的 新 雪 雪 有 衣 。 调 查 时 研究 区 内 雪 朋 已 发 生 且 仍 处 
于 灾害 孕育 阶段 , 雪 层 裂缝 有 加 剧变 形 趋势 ,风力 
作用 下 雪 榴 自重 增 大 ,有 超过 雪 的 抗 断 强度 的 趋 
势 ,整体 稳定 性 差 , 在 降雪 、 风 等 外 力作 用 下 具有 发 
生 大 规模 雪 衣 的 可 能 性 。 

(3) 以 坡 面 上 方 积 雪 平 台 为 潜在 释放 区 的 坡 面 
型 雪 朋 ,释放 量 可 达 8.2669x10' m’, 运 动 时 长 约 为 
128 s, 并 在 120 s 内 堆积 区 流动 高 度 达 到 最 大 , 约 为 
3.55 m, 流 动 速度 最 快 可 达 18.34 ms ,冲击 力 达 到 
32.67 kPa。 而 堆积 量 达 到 1.8525x10:m, 影 响 面 积 
在 3369.7 mx。 通 过 相互 验证 , 坡 面 型 雪崩 并 非 积 雪 
平台 的 释放 ,地 面 调查 结果 与 数值 模拟 结果 存在 
差异 。 

(4) 沟 模 一 坡 面 复合 型 雪 朋 在 不 同 断 裂 深 度 下 
进行 模拟 ,其 堆积 体形 态 受 坡 面 影响 均 形 成 中 间 
高 .四周 低 的 连续 性 椎 状 堆积 体 ,与 现场 地 面 调查 
和 航 摄 解 译 堆 积 形 态 结果 一 致 ; 而 在 雪 层 断裂 深度 
减少 40% ,平均 堆积 深度 与 现场 堆积 体 测量 结果 相 
近 ; 因 此 , 沟 槽 一 坡 面 复合 型 雪 朋 为 沟 槽 坡 面 雪 层 
断裂 释放 旦 断裂 深度 仅 为 临界 厚度 值 的 60% 左 
右 。 雪 层 断 裂 深 度 减 少 40% 下 雪崩 持续 时 间接 近 
于 300 s, 并 形成 1178.5 mm 堆积 区 ,平均 堆积 深度 为 
1.64 mm, 堆积 量 3107.76 n» ,堆积 区 最 大 流动 速度 可 
356.58 ms ,最 大 冲击 力 为 17.97 kPa。 
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Numerical simulation of avalanche process in Aerxiangou, West Tianshan 


Mountains, based on air-ground cooperative investigation 


ZHANG Tianyi", LIU Jie", YANG Zhiwei’, WANG Bin’, CHENG Qiulian™ 
(1. School of Transportation and Logistics Engineering, Xinjiang Agricultural University, Urumqi 830052, Xinjiang, 
China; 2. Xinjiang Transportation Planning, Survey and Design Institute Co., Ltd., Technology Research and 


Development Center, Urumqi 830006, Xinjiang, China) 


Abstract: This study was implemented to accurately identify the avalanche flow characteristics and flow 
information, and comprehensively analyze avalanche motion. This study was based on UAV tilt photography 
technology to obtain high-resolution aerial photography data, taking the avalanche-prone area of Aerxiangou as 
an example. Through on-site investigation and UAV remote sensing interpretation to detect avalanche activity in 
a high-resolution manner, the goals were to determine the input parameters of the RAMMS model, to simulate 
and reconstruct different types of avalanche events on this basis, and to comparatively analyze the differences 
among the results of conventional ground-based investigation, UAV remote sensing interpretation, and simulation 
results to explore avalanche activity in different types and different snow layer release conditions. The results of 
the study show that (1) the avalanche investigation and analysis system based around tilt photography technology, 
which combines conventional ground- based investigation methods with UAV remote sensing and numerical 
simulation to verify each other, improves the accuracy of the assessment of disaster development status. (2) In 
mid- February, the snow on the slopes of Aerxiangou approaches the critical thickness value, and continuous 
snowfall destabilizes the snow layer and triggers new avalanches. The investigation is still in the disaster 
breeding stage, the snow layer cracks intensified deformation, the role of the wind snow eave self- weight 
gradually increased, there is more than the trend of the breaking strength of the snow, and the overall stability is 
poor. (3) In slope-type avalanches with a snow platform above the slope surface as the potential release area, the 
release volume can reach 8.2669 x 10* m°, the movement duration is about 128 s, and the flow height of the 
accumulation area peaks in 120 s at about 3.55 m, the flow velocity is about 18.34 m.s ', and the impact force is 
about 32.67 kPa. In addition, the accumulation area is formed into an accumulation with an area of 3369.7 m and 
a volume of 1.8525x 10 m! of the pile. Through mutual verification, the slope-type avalanche does not involve 
release of the snow platform, and there is a discrepancy between the ground- based investigation results and 
numerical simulation interpretation results. (4) Trench-slope composite avalanches are released by fracture of the 
snow layer on the trench-slope, where the depth of fracture is only about 6096 of the critical thickness value, the 
avalanche duration is close to 300 s in this case, and the impact range is 1178.5 m in the accumulation area, with 
an average accumulation depth of 1.64 m. The flushing-out volume is 3107.76 m’, the maximum flow rate in the 
accumulation area is 6.58 m-s ', and the maximum impact force is 17.97 kPa. The results of the ground-based 
investigation are roughly the same as those of the numerical simulation based on the 3D model. The results of the 
study have improved the accuracy with which avalanche event information can be acquired and can provide 
strong data support and a scientific basis for predicting future avalanche potential hazards, risk avoidance, and 
disaster emergency response. 
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